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Abstract 

Current-driven vortex wall dynamics is studied by means of a 2-d analytical model and micro- 
magnetic simulation. By constructing a trial function for the vortex wall in the magnetic wire, we 
analytically solve for domain wall velocity and deformation in the presence of the current-induced 
spin torque. A critical current for the domain wall transformation from the vortex wall to the 
transverse wall is calculated. A comparison between the field- and current-driven wall dynamics is 
carried out. Micromagnetic simulations are performed to verify our analytical results. 



PACS numbers: 75.60.Ch, 75.75.+a, 75.70.Kw 



I. INTRODUCTION 



Magnetic domain walls in magnetic films have various structural forms which are de- 
termined by geometrical and material parameters. For a magnetic wire with an "infinite" 
length and a finite width, there are commonly two types of domain walls: transverse wall 
(TW) and vortex wall (VW). Both of them are stable. Depending on the wire thickness 
and width, one of the walls is usually more stable . However, in a certain range of the 
parameters, the static energies of these two walls are comparable and thus one can produce 



both types of wal 
create either wall 



s in the same wire . By using different initialization methods, one can 
l|. When a magnetic field or an electric current is applied, both TW and 
VW are able to move along the wire. The dynamics of the walls is generally very complex 
and micromagnetic simulations are required in order to describe the details of the domain 
wall motion. For the TW, a simplified and yet very insightful analytical treatment was 
developed by Walker Q|. A 1-d wall profile, i.e., the magnetization direction in the wall 
depends only on the coordinate along the wire, has been used to approximate the TW pro- 
file. With this approach, one can analytically calculate the wall velocity and wall distortion 
in the presence of magnetic field and electrical current ^J. For the VW, however, the 1-d 
wall profile fails to capture the wall structure and one needs at least to use a 2-d model to 
approximately characterize the vortex structure. In this paper, we propose such a 2-d model 
for the VW. Our focus will be on the analytical calculation of the dynamic behavior of the 
VW. Within our model, we are able to describe the vortex wall motion, including the wall 
distortion, wall velocity and wall structure transformation, in terms of material parameters 
and external magnetic field or electric current. In particularly, we show how the vortex core 
moves toward the edge of the wire when a current or a field is applied along the wire. With 
a sufficiently large current or field, the vortex core may vanish at the boundaries of the wire 
edges and the transformation from the VW to the TW occurs. This paper is organized as 
follows. In Sec. II, we develop the analytical 2-d model for the VW. An equation of motion 
for the domain wall is established. The steady state motions driven by the current and by 
the field are investigated. We also compare the dynamics between the TW and the VW. 
In Sec. Ill, the micromagnetic simulations are performed. We compare the simulated results 
with the analytical ones. Finally, we summarize the different features of dynamics for TW 
and VW in Sec. IV. 
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II. ANALYTICAL MODEL 



A. Equation of motion 

We start with the generalized Landau-Lifshitz-Gilbert equation including the spin transfer 
torque terms 0]: 

— — = -7M x H eff H M x — — (1) 

^ ' eff M s dt K J 

-M x M x — — -M x 



M? V dx J M s dx 

where 7 is the gyromagnetic ratio, and H e ff — ~wi * s ^ ne e ^ ec ^ ve magnetic field, and the W 
is the total energy density which could be written as W = A'V 2 m+(K / Af/)(M x e x ) 2 — H e - 
M— (l/2)Hd-M, where ^4 is exchange constant, K is anisotropy, H e is external field, and Hd 
is magnetostatic field, a is the Gilbert damping parameter, and bj = Pj e f^B / eM s (l+^ 2 ) and 
cj = £bj, where P is the spin polarization of the current; j e is the current density along the 
length direction of magnetic wire; \xb is Bohr magneton, M s is saturation magnetization and 
£ (small, ~ a) is a dimensionless constant which describes the degree of the nonadiabaticity 
between the spin of the nonequilibrium conduction electrons and local magnetization. 

To describe the motion of an entire domain wall, it is useful to introduce a total force 
acting on the wall. Following Thiele p, we define the total force 



F{9,(f)} = J dVVW = J 



8W 5W 
(W) + ttW) 



dV (2) 



59 y ' 8(f) 

where 9, <fr are the angular components of M in the spherical coordinate. For the steady- 
state motion of a domain wall, we may write 9 = 9{j — vt),(j) — 0(r — vt), where v is the 
steady velocity, then we have 

e = - V -V9, j> = -vV<l>. (3) 

The above steady-state condition immediately reduces the temporal-spatial differential equa- 
tion, Eq. (1), to a differential equation with spatial variables only. By writing Eq. (1) in the 
angular components and by placing them into Eq. (2), we obtain the equation of motion for 
the domain wall Q|: 

F + G x (v+ 6jx) + V ■ (at; + cjx) = (4) 
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where G is the domain gyrocoupling vector 

G = -M^- 1 J dVsm6{V6 x V0) (5) 

and V is the domain dissipation dyadic (tensor) 

V = -M s j~ 1 J dV(X79X79 + sin 2 9X7 <pV(f)) (6) 

The domain wall force defined in Eq. (2) can be simplified when the domain wall undergoes 
uniform motion. Let us separate the force in terms of the internal force F m and the external 
force F ext , 

F = F in + F ext . (7) 

F m contains all the forces from the internal energy including anisotropy energy, exchange 
energy, and the magnetostatic energy. When one sums over the internal energy contribution, 
the total internal force vanishes due to Newton's third law. Therefore, one may simply 
consider the external energy contribution to the force in Eq. (4) and (2), i.e., 



text 



dV 



H M). (8) 



where we have assumed that the external energy is solely from the external field, W = 
— H • M. We will show later that we must consider other external forces on the vortex wall 
when the wall reaches the boundary of the wire. Note that if the profile of domain structures 
(i.e. M(x, y, z) in the moving frame of the steady motion) is determined, the gyrocoupling 
vector, dissipation dyadic and static force can be calculated from Eqs. (5), (6) and (8), and 
the steady velocity will be then readily derived from Eq. (4). 

In this Section, we shall apply the equation of motion, Eq. (4) , to study the domain wall 
velocity of the vortex wall, driven by an external magnetic field and by spin transfer torques. 
Let us first consider a simplified head to head transverse wall as shown in Fig. (la). For the 
transverse wall, we assume the wall profile can be modeled by the Walker's trial function 

Q 

<f>{x) = 2tan" 1 exp , 9{x) = ~. (9) 

where <p(x) and 9(x) are the angles between the direction of the magnetization and the wire 
length direction (+x-axis) and wire plane normal (+z-axis) respectively. A is the domain 
wall width. By placing the above wall profile into Eq. (5)- (8), we find that the gyrocou- 
pling vector is zero G = and the dissipation dyadic has only one non-zero component 



T^xx — —2M s / , yA. The external force is F = 2HM S £. By inserting them into Eq. (4), we 
immediately find the velocity: 

* = 2™-£i. (io) 



This result had been obtained previously 

It is difficult to analytically model the profile of the vortex wall depicted in Fig. (lb) by a 
single elementary function as we did for the transverse wall. From previous simulation results 
Q,|2|, the VW structure might be understood as two symmetrical transverse walls diagonally 
crossing the wire and a central vortex core connecting the two TWs. It is noted that the two 
TWs have opposite polarities, namely, the magnetization of the centers of these two TWs 
orient in opposite directions. The transitional region between these two TWs is sometimes 
called as the Bloch line [Uj , which characterizes the wall-polarity reversal in analogy to the 
Bloch wall. For the VW we consider here, this transitional region contains a vortex core. The 
magnetization of the inner vortex core has a significant out-of-plane component and thus 
the inner radius of the vortex core must be very small since the out-of-plane magnetization 
enhances demagnetization energy ^| . Outside the inner vortex core the magnetization lies 
in the plane and the outer radius of vortex core is limited by the transverse wall width and 
the wire width. To characterize the entire VW profile, we separate the wall into three parts: 
two transverse walls and a vortex core (schematically shown in Fig. (lc) and (le)), and for 
e) : they will be assigned to different trial functions given below . For 



the model in Fig, 
the vortex part [12, 



ej, tneywi. 



2 tan" 1 



2 - 



, (0<x 2 + y 2 <r 2 core ), 
(r 2 core <x 2 + y 2 <R 2 ), 

,2 



(11) 



(j) = q ■ arg(x + iy) + c|, (0 < x 2 + y 2 < R 2 ). 

where r core and R are the inner-core and outer radius of the vortex respectively, q(= 
±1, ±2...) is the vorticity of the vortex, c(= ±1) is the chirality of the vortex, and % = a/— T. 
Here we just use the arg function of complex variable x + iy as a convenient way to express 
[3]. In this paper, we only consider a single vortex (q = 1). And we note that the 
vortex profile we introduced in Eq. (11) does not include its image profile Q|. For the wire 
structure, the image profile would consists of a series of terms formed by multiple reflections 



of two wire boundary planes. Including these images would make the analytical calculation 



intractable. Instead, we introduce a 'restoring force' |18J] due to the induced charges along 
the boundaries of the nanowire. This restoring force has similar roles of the 'image vortices' 
by the method of image [15]. For the two transverse walls ^,0] where (—00 < x < +00), 



9 = tt/2, 



-yo<y< yo), 



2tan- 1 exp (2^0) , (-y -5y<y<0), 



(12) 



2 tan- 1 exp (^f ) , (0 < y < y - Sy) 

where Xq is the distance from the centers of each TW part to the vortex core center, w' 
and w" are the wall widths of two TWs, yo is the half width of the magnetic wire, 5y is the 
displacement of vortex core away form the wire center in the y-direction. We realize that 
the TW and vortex core defined above are not strictly valid because the TW is extended to 
the vortex region when < x 2 + y 2 < R 2 . However, this ill-defined overlapped region does 
not contribute significant errors since V6 1 and is exponentially small in this region. By 
using the above wall profile, we are able to calculate G, V and F defined in Eqs. (5), (6), 
and (8) in a similar way as the calculations in j^, j|. The steady-state velocity can then be 
calculated from Eq. (4). 



1. Current- driven domain- wall motion 



We consider the current-driven vortex wall motion without the external field so that 
F = 0, see Eq. (8). Let us apply the equation of motion, Eq. (4), to the vortex part of the 
wall and neglect its interaction with two TW parts. By placing Eq. (11) into Eq. (5) and 
(6), we find the explicit expressions for G and V: 

G = -lirMs^pi = G v z (13) 

R 

V = -M s 7 _1 7rln— (xx + yy) (14) 

(X 

= P„(xx + yy) 

where x, y, z are unit vectors in the directions of the wire length, width and thickness, 
respectively. The p(— ±1) represents the polarity of the vortex core (#(0) = 0, or, tc), and 
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a is the lattice constant of the crystal structure of wire material. Note that in the calculation 
of T>, we have neglected a small part of the vortex core (~ a 2 ), and this does not considerably 



effect the wall dynamics [8|. By placing Eqs. (13)-(14) and F = into Eq. (4), we have 



-G v v v + VJav x + cj) 



G v {v x + bj) + aV 

v v y 



(15) 



and one immediately finds: 



XQ 



(16) 



It is shown that the velocity of VW has two components: v xg along the wire and v yo per- 
pendicular to the wire. The nonzero (except £ = a) is caused by the gyrotropic term 
(nonzero G) of the vortex dynamics 6]. The sign of v yo (to which edge of the wire) is 
determined by the sign of the product p(a — £), i.e., it is vortex polarity p dependent. If we 
characterize the vortex by a Bloch line, the perpendicular motion might be understood in 
terms of the moving Bloch line inside and along the wall surface jl 1| . 

We recall that the velocity of the TW is v x = —cj/a and v y = in the steady-state 
motion. The different velocities for the vortex core and the transverse wall can not be a 
steady state solution of the entire wall. In fact, the interaction between these different parts 
of the wall must be included in order to reach a common velocity in the steady state. One 
can immediately see that the motion of the vortex core along the ^-direction would push 
one TW and pull the other. As a result, one TW expands and the other gets compressed, 
see Fig. 1(e). Conversely, these distorted TWs would produce a reacting force to the vortex 
core; sometimes it is called a restoring force QQQ. If we model the forces between the 
vortex core and the TW by elastic potentials in analogy to spring-connected particles, we 
must added the spring force as an external force when one applies Eq. (4) to each individual 
parts of the vortex wall. Then, two scenarios are possible: the reacting force is strong enough 
to completely halt the perpendicular wall velocity; or the reacting force is unable to stop the 
vortex from colliding to the wall edge. In the former case, the perpendicular wall velocity 
eventually reaches zero and a steady-state wall velocity along the wire is achieved. The final 
wall velocity is precisely same as that of the transverse wall (i.e., v x = —cj/a and v y = 0, 



see calculations below). In the latter case, the vortex core collides with the wire edge and 
vortex core can either vanishes (move out of the wire) or be reflected, i.e., the transformation 
of the vortex wall to other types of the walls occurs. Usually, the vortex core disappears at 
the edge and the wall becomes a TW 

When a steady-state motion is achieved, we may apply Eq. (4) to the entire wall. The 
interaction between the vortex core and the TW walls becomes internal force and it does 
not enter Eq. (4). However,there is an external force acting on the entire wall, which is due 
to the magnetic charges coming from the distorted TW parts at the boundaries. This force 
can be again analogy to the counteractive force from the fixed boundaries to the compressed 
springs. If we model this force linearly depending on the displacement of the vortex core Sy 
in the y-direction, i.e., F re = —n5y where the constant k is discussed in Ref. [l^ , we find, 
from Eq. (4), 

-G v v v + V x Jav x + cj) = 



(17) 



-K5y + G v (v x + bj) + aV yy v y = 
where we have used and determined 



G = G v z (18) 

V = (p v + V t )xx + V v yy (19) 
= T> xx xk + Vyyyy 

where T> t = —2M s , -f~ 1 ( yo ~f y + yo *f y ) . Note that the TW parts have no contribution to the 
gyrocoupling vector G, however they modify the dissipation dyadic with T> t . 
In the steady state of motion, v y = 0, we immediately find from Eq. (17) 

v x = (20) 



5y = ^(l-^) (21) 



and 

8y — — 

k \ a t 

The above simplified analysis shows that the steady state domain wall velocity is a uni- 
versal constant v x = —cj/a, independent of the wall structure, as long as the restoring force 
from the edges of the wire is strong enough to halt the wall motion in the direction perpen- 
dicular to the wire. We can qualitatively estimate the critical current 6} when the restoring 
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force just barely prevents the wall colliding with the wire edge, i.e., when by = yo (yo is the 
half width of the wire). By using Eq. (21), we have 

G v \ a 

If the current is larger than the above current density, the restoring force is unable to 
stop the vortex core colliding to the boundary and the VW transforms into a TW. The 
wall transformation may be qualitatively understood in terms of the conservation of the 
topological charges 16]: when the vortex core reaches the edge, the vortex with the winding 
number +1 is absorbed by one of the edge defects (one of the transverse wall parts) with 
winding number —1/2, then a +1/2 edge defect must be created to conserve the topological 
charges. This +1/2 edge defect together with the —1/2 edge defect on the other side of wire 
edge is equivalent to a simple transverse wall with a uniform polarity. 




2. Field-driven domain-wall motion 



By adding the external force due to the magnetic field (defined in Eq. (8)) to Eq. (17), 
we can similarly derive the domain-wall steady motion velocity. In the steady state, the 
external force on the entire wall is 



Ay HM s yi + cnHM s (R - a)y = F x x + F y y. 



(23) 



In the absence of the current where cj = bj = 0, we find that from Eq. (17) the steady-state 
velocity along the wire (v y = 0) is 



lHy 



a 



4 



n l n R + 2(yo^y + y»+ty) 

a V in' hi" I 



(24) 



and the displacement of the vortex center in the direction of the wire width is 

8y = y - 



(25) 



As in the case of TW, the wall velocity is not universal and it depends on the details of 
the wall structure. Comparing with TW, the velocity of the VW is smaller n|; this can 
be qualitatively understood in terms of a smaller wall width for the VW. Let us take the 
width of the TW part Q] as vo ~f y = m+ ^ v = 7r (again, yo is half width of the wire) . And 
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the width of the TW (A in Eq. (9)) in the same wire is also taken as A = 2y^j-K. For the 
VW, we may define an effective wall width from Eq. (24) 

«W = A (^T^J (26) 

For example, for a wire R = w/2 = 20nm and the lattice constant a = 0.5nm, the ratio of 
the TW to VW velocity is about 2 + \ In f = 3.8. 

When the magnetic field is large, the restoring force is unable to stop the wall colliding 
with the wall edge. The vortex wall will be annihilated and the transformation to the TW 
occurs. Similarly, we may estimate the critical field H c e for the wall transformation by the 
condition by = yo- We thus find 

(27) 

III. MICROMAGNETIC SIMULATION 



«3/0 



lG v w v 



a 



+ ckM s {R-cl) 



A. Simulation procedure 

To verify the analytical results derived above, we perform micromagnetic simulations by 
directly solving Eq. (1) for a defect-free magnetic wire sample. The geometrical size of the 
wire is 2/xm long (x direction), 128nm wide (y direction) and 8nm thick (z direction). The 
grid size is taken as 4 x 4 x 8nm 3 . The material parameters are: the exchange constant 
A = 1.3 x 10~ 6 erg/cm, the anisotropy field = 0(Oe), the saturation magnetization 
M s = 800emu/cc, the spin polarization P = 0.5, the damping parameter a and the non- 
adiabaticity factor £ varying from 0.01 to 0.08. 

The magnetization at both ends is set to be along the x-direction and direct inward to 
the wire and we use free boundary conditions for other dimensions. Since the domain wall 
may move several micrometers in some cases, it is important to keep the domain wall far 
away from the ends to reduce the influence of magnetostatic fields from the wire ends: we 
shift the x-coordinate of the entire domain wall after each numerical iteration so that the 
center of the wall is always located at the center of the wire. By using two initial domain 
configurations similar to pj, we generate two types of walls in the middle of the wire, the 
TW in Fig. (la) and the VW. The VW has four equivalent configurations corresponding to 
different polarity (p = ±1) and chirality (c = ±1). We only show the VW for p = —1 and 
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c = 1 in Fig. (lb). The TW and VW shown in Fig. (1) are our initial wall configurations. 
At time t — 0, a constant spin torque or a magnetic field is applied to the wire, and we 
calculate the wall motion as a function of time until the wall reaches steady motion or until 
a sufficiently long time that the wall motion remains non-uniform. 



B. Simulation results 



To compare the analytical results on the universal domain wall velocity, we first consider 
applying a sufficiently small spin torque (bj = 25(m/s)) so that the vortex core remains 
inside the wire. We vary damping parameter a and non-adiabaticity coefficient £ from 0.01- 
0.08 to determine the relation between the domain wall velocity and these two parameters. 
In Fig. (2), We have found that the initial velocity of the VW in the x-direction is always 
bj (independent a and £), while the total displacement of the center of the wall in the 
y-direction quantitatively agrees with Eq. (21). As expected, the ^-component velocity 
gradually reduces to zero when the center of the core moves toward the edge. The restoring 
force eventually stops the wall motion in the ^-direction, and a steady state motion with 
the velocity v x = —cj/a along the x-direction is reached. In Fig. (3), we show the relations 
of the terminal velocity v x with both a and £. These simulated results are in excellent 
agreement with our analytical formula. 

When we increase the applied current density (or the strength of the spin torques) above 
a critical value, the simulations show the vortex core moves outside of the nanowire from 
;he edge and the VW transforms to TW. The transformed TW moves at the velocity —cj/a 
2o| . We note that the critical current density depends on the parameters. A special case 
is when £ = a. We find that the VW can maintain its wall structure up to a much higher 
current density. This is because the vortex wall does not move toward the wire edge, see 
Eq. (21). However, this accidental case should not be considered as a general property of the 
VW motion. When the current density is further increased, the TW will also deform during 
its motion and the velocity is highly non- uniform, known as Walker's breakdown 3j. The 
alternative appearance of transverse wall and vortex wall (or anti- vortex wall) |l9] is seen 
in simulations (not shown here). In our analytical calculations, we exclude these extremely 
complicated wall dynamics. 

As a comparison, we also simulate the field-driven domain wall motion at sufficient small 
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fields. The behavior of the VW driven by field is similar to the current- driven case: namely, 
besides the motion along the wire, the vortex core moves toward the edge and eventually 
stops due to the restoring force. At a critical field, the VW overcomes the restoring force 
and the VW transforms into TW. While the steady velocity of TW and VW are both 
proportional to the magnetic field H and a -1 , see Fig. (4), the VW has a smaller velocity 
compared to the TW for the same field. 

Finally, we point out that the analytical and simulation analysis in this paper is for 
a defect-free wire. In a realistic wire, the defects are unavoidable. We expect that our 
analytical results may not be applicable. For example, the steady state wall velocity may 
not be universal and the defect pinning dep ends on the detail of the wall structure. Indeed, 
we have seen in our earlier simulation [20| that the critical current is larger for depinning 
a TW than for a VW. An even more remarkable finding is that the VW firstly begins 
depinning, then transforms to a TW, which is finally pinned by defects again (wall stops). 
The phenomenon occurs when one applies a moderate current density between the critical 

n 

values of two depinning currents of the VW and TW 20]. The above feature has been 
captured in recent experiment jlTI ] . 



IV. CONCLUSION 



In summary, we have calculated the velocity of the VW in the magnetic wire. In a 
defect-free wire, the terminal velocity is independent of the wall structure for the current- 
driven steady motion, while for the field-driven case, the effective domain-wall width is 
a key parameter for the domain-wall mobility. The transformation between the VW and 
the TW is explained as the consequence of the perpendicular motion of the vortex core 
and the conservation of topological charges. Our model is further supported by numerical 
solutions and in agreement with experimental results. The research was supported by NSF 
(DMR-0314456 and ECS 0223568). 
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Figure Caption 



FIG.l Magnetization patterns of (a) a transverse wall and (b) a vortex wall, without an 
applied current or field, (c) Schematic model for the vortex wall (b): a vortex wall is 
modeled by two TWs and a vortex core, shown in gray areas. The vortex core is at the 
center of the wire and ±x are the positions of centers of the TWs. R is the outer radius of 
the vortex core, w is the wall width of two TW walls, and yo is the half width of the wire, 
(d) The magnetization pattern of the vortex wall in the presence of the current, which is 
calculated by micromagnetic simulation in the Sec. III. The parameters are: the spin torque 
bj = 25(m/s), and damping parameter a = 0.05, non-adiabaticity coefficient £ = 0.04. (e) 
The schematic model for the vortex wall (d). 5y is the displacement of the vortex toward 
the edge of the wire. The wall structure is not symmetrical, and one TW is expanded and 
the other is compressed. 

FIG. 2 The dependence of (a) x-component of initial velocity v XQ and (b) the y-component 
displacement 8y as a function of the nonadiabaticity parameter £ for two different polarity 
VWs (p = ±1). We have used bj = 25(m/s) and a = 0.05. The straight lines are simple 
data fits. 

FIG. 3 Steady-state velocity of TW and VW as a function of (a) the damping parameter, 
and (b) the non-adiabaticity parameter. The current density is bj = 25m/ s. The straight 
lines are simple data fits. 

FIG. 4 Steady-state velocities of both TW and VW as a function of (a) the magnetic field 
and (b) the damping parameter. The current density is bj = 25m/ s. The straight lines are 
simple data fits. 
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